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Abstract 

Demanding the analyticity of hadronic observables (calculated in terms of power series of the 
running coupling) as a whole, we show that they are free of the Landau singularity. Employing 
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to two different hard-scattering functions in perturbative QCD: the electromagnetic pion form 
factor to leading order and the inclusive cross section of the Drell-Yan process. In the second case, 
the leading nonperturbative power correction in 6Aqcd gives rise to a Sudakov-like exponential 
factor in the impact parameter space which provides enhancement rather than suppression. 
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1. Introductory remarks 

One of the most crucial questions in comparing QCD with experimental data is whether 
predictions derived at the parton level can claim validity also at the hadron level. Would we 
be able to calculate binding (confinement) effects of quarks and gluons reliably, this prob- 
lem would be only marginal because then calculated (i.e., partonic) and physical observables 
would almost coincide. Lacking such a theoretical scheme, the strategy must be to make per- 
turbative predictions less sensitive to nonperturbative effects. Since parton binding effects 
become important at large distances (i.e., small momenta), we must ensure that sensitivity 
of hadronic observables to the large- distance domain is minimized (infrared safety). For 
such infrared-improved observables, one may expect that deviations between theory and ex- 
periment, originating from the differences between partons and hadrons, become practically 
irrelevant. 

Normally, large-distance effects are factorized into universal, wave-function parts that 
cannot be treated perturbatively. However, infrared (IR) sensitivity may also reside in the 
initially separated hard part that should describe partonic subprocesses involving by defi- 
nition hard propagators and hard vertices. In fact, beyond leading order, this assumption 
cannot a priori be satisfied because the momentum flow in such Feynman diagrams may 
become much smaller than the large scale of the process, say, the momentum transfer Q 2 . 
More precisely, the average gluon momentum < k > flowing inside a partonic subprocess 
decreases with the order of the perturbation expansion and may eventually become pro- 
portional to Aqcd, implying that the resummed perturbation series is defined only up to a 
certain power accuracy. To make the result of the perturbative calculation unambiguous, 
one has to compensate this power correction in the perturbative sector by a corresponding 
term of the same form, originating from the nonperturbative regime. 

In processes which involve the emission of virtual gluon quanta of low momentum, the 
strategy must be to resum their contributions to all orders of the strong coupling constant. 
This gives rise to exponentially suppressing factors of the reaction amplitude (or cross- 
section) of the Sudakov type with exponents containing double and single logarithms of the 
large mass scale. However, because of the Landau singularity of the running coupling at 
transverse distances b oc 1/Aq C d (where b is the impact parameter conjugate to the parton's 
transverse momentum fcjj, an essential singularity appears in the Sudakov factor. Thus, 
one has to consider power corrections of O (6 2 Aq CD ) , which, though negligible for small b 
relative to logarithmic corrections In (6 2 Aq CD ) , may become important for larger values of 
the impact parameter. 

In this letter we will describe a general methodology to treat power series in the running 
strong coupling. To be more precise, we will address these questions having recourse to 
two processes: one to which the OPE applies, namely the pion electromagnetic form factor 
at leading perturbative order, and another, the Drell-Yan process, to which the OPE is 
not applicable. The first one is a typical example of an exclusive process with registered 
hadrons in the initial and final states (for a recent review, see, e.g., pj). The Drell-Yan 
mechanism, on the other hand, has two identified hadrons in the initial state and a lepton 
pair (plus unspecified particles) in the final state, whose transverse momentum distribution 
is proportional to the large invariant mass of the materialized photon. 

Our goal in the second case is to obtain not only the usual resummed (Sudakov) expres- 
sion, comprising logarithmic corrections due to soft-gluon radiation, but also to include the 
leading power correction as well, specifying, in particular, its concomitant coefficient. This 
becomes possible within a theoretical scheme, which models the IR behavior of the run- 
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ning coupling by demanding analyticity of physical observables (in the complex Q 2 plane) 
as a whole - as opposed to imposing analyticity of individual powers, i.e., order by order 
in perturbation theory -, while preserving renormalization-erroup invariance (references and 
additional information can be found in the recent surveys [2, H[ ) . 

Our own physical viewpoint can be summarized as follows: One option to access non- 
perturbative effects is to change degrees of freedom and replace QCD by some low-energy 
effective theory. The other option, and that actually adopted here, is to retain the usual 
QCD degrees of freedom, but to demand that hadronic observables, calculated perturbatively 
with them, are analytic as a whole in the Q 2 plane. This analytization condition entails a 
singularity-free expression for the strong running coupling in both the spacelike and the 
timelike region, rendering hadronic observables IR safe. In this way, we are able to calculate 
explicitly in our second example the 6 2 Aq CD power correction to the Drell-Yan cross-section 
and show that after exponentiation it amounts to a Sudakov type factor that can compete 
in magnitude with the resummed double logarithms due to gluonic radiative corrections be- 
cause it has the reverse sign, whence providing enhancement rather than suppression. The 
properties of this new factor are addressed ^and its (almost) Gaussian dependence on b is 
discussed and compared to previous works 



2. Power corrections to pion form factor 

Following the rationale of analyticity in the sense just described, the leading-order fac- 
torized pion form factor reads ja [7| 



[Q 2 F n (Q 2 )] an = JdxJ dy [$° ut (y,Q 2 ) T H (x,y,Q 2 ,a s (q 2 )) ^ {x,Q 2 )]^ , (1) 

where (x, Q 2 ) is the process-independent pion distribution amplitude, encoding the non- 
perturbative binding dynamics of the valence quarks at the resolution scale Q 2 , each carrying 
light-cone momentum fractions x\ — x (quark) and 22 = 1 —x = x (antiquark) (xi = kf / ' P + ) 
of the pion's momentum P^, and the superscripts in and out denote, respectively, incident 
(incoming) and final (outgoing) intact pions. 

Employing the asymptotic pion distribution amplitude, evolution effects can be neglected 
(H and the analyticity requirement resides only in the hard-scattering part. Hence, we have 

[Q 2 F n (Q 2 )] an = A^ dx J dyxyxy [t h (x,y,Q 2 ,a s (Q 2 ))}^ , (2) 
where to leading perturbative order, the hard-scattering amplitude is given by 

'2a s (Q 2 xy) la s (Q 2 xy) 



3 xy 3 xy 



(3) 



and the constant A takes account of the correct normalization of the pion distribution 
amplitude. Had we not imposed the requirement of analytization, expression (J2J) would lead 
to an asymptotic series in the coupling constant a s (Q 2 ) that is not Borel summable, as 
noticed by Agaev (see also [Io|). 

Global analytization 0, El (see also jlil . ) now means that 
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where the spectral density p( n '(£) is the dispersive conjugate of all powers n of a*. For the 
leading-order expression under consideration the spectral density is 



p (Q 2 ) = lma s (-Q 2 ) 
with A = Aqcd- Then Eq. (J3j) reduces to 

M<? 2 )L = ^/ c 



7T 



A In 2 g + 7T 2 



(5) 



2; - Q 2 + ie 



a s {z) 



(6) 



where C is a closed contour in the complex z-plane with a branch cut along the negative 
real axis, so that = Po/Air with /3 = 11 — 2Nf/3) 



Me 2 )L 



1 / 1 



Pi \ In 



1 - 



(7) 



an expression recently proposed by Shirkov and Solovtsov [1 
Recasting now the strong coupling in the form 



a s (z) 



1 1 



Pi In 



A 2 



± / dcrexp (=F<j/?i lnA 2 /^) 



(8) 



with the plus sign corresponding to the case |-z|/A 2 > 1 and the minus one to |^|/A 2 < 1, 
and inserting it into Eq. (J2J), we find after some standard manipulations the Borel transform 
of the scaled pion form factor at leading perturbative order: 



[q 2 f w (qr 



2M(1) 



daexp {-a(3i \nQ 2 / K 2 )tx(g 



,(i) 

an 



(9) 



Here the Borel image of the form factor reads 



7r (a 



l6irC F A 



sin (7ij3i(j) 



71 



A 2 
Q 2 



dx / dyxy 
Jo 

2\ 2o-/3i 



Q_ 

A 2 



A 2 
Q 



(10) 



We stress that this expression has no IR renormalons in contrast to approaches which use the 
conventional one-loop a s parameterization (see, e.g., |H,[l3]). Carrying out the integrations, 
we then obtain the following final result 



7T 



a){]} = 16nC F A 



sin 



(irfac) /A 2 V^ 1CT f Q2/A2 dw 



x 



2 F X l,l;l + Aa; 



w 



+ 



l + w T V. Q 



2 F X 1,1;2 + ^; 



1 + w 



1 + wJ 1 + Act 

with 2F1 being the hypergeometric function and 4>(w) denoting the abbreviation 
4>(w) = —(1 + w) \nw — 2(1 — w) (<K^) > when w < 1) 



;i2) 
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Inserting this expression into Eq. Q, the integration over the Borel parameter a can be 
performed without any ambiguity to arrive at the following result for the pion form factor 
(in leading order) 



[Q 2 F W (Q 2 )]^ = IQttCfA^- [ dwcf>(w) 

Pi Jo 



A 2 



(13) 



Hence, the remaining integration can be carried out analytically and the final result is 



[Q 2 f„(Q 2 )] 



2M(1) 
an 



13 A 2 



lnlgl+2 



A 4 

A 2 / A 2 

+ q 2 ~{ + q 2 ~ 



A 2 " 
A^ 



T I ^ 

A 2 



Ul ' A 2 



2 ln A 2 " 



A 2 A^_A 2 _A A"_ 

" Q 2 + Q i Q 2 V + Q 2 ~ 



n=l 



1 /A 2 



2 \ n 



In 



1 



7> z 



1 - 



A2 \ "1 

Q 2 



A 2 



(14) 



It is important to notice that the above expression remains analytic all the way down to the 
limit — > 1. This is to be contrasted with the corresponding result found by Agaev in 
(his equation (16)) which comes out divergent in this limit and has to be regularized. 
For the physically relevant case Q 2 3> A 2 , Eq. (fTl|) simplifies to 



\a s (Q 2 ) + O (a 2 ) 



1 A 2 r 



71 



+ 3 



(15) 



where LI2 is the di-logarithm (or Spence) function, defined by 

" x dt 
n hit 



Li 2 (x) = — 



(16) 



The computation carried out above can be extended to any desired order of the hard- 
scattering amplitude Tr. The resulting expressions are always well-defined without the need 
of employing any (additional) IR regularization (3). 



3. Power corrections to Drell-Yan process 

As a second example of our methodology, we discuss the derivation of power corrections 



5 



to the inclusive Drell-Yan cross-section with the large scale Q 2 being here the invariant 
lepton pair mass. Power corrections to this process with the help of renormalons have been 
discussed in jHESEilE, 13 • Let us start our discussion here by evaluating the logarithmic 
derivative 1 of the unrenormalized expression for the eikonalized Drell-Yan cross-section |23| 
(see also 24]|) after exponentiation of infrared divergences to the lowest order of perturbation 
theory, adopting here and below for the ease of comparison the notations of Korchemsky 
and Sterman 



d In Wby 
dlnQ 2 



AC F fx 



2c 



d 2 ~ 2e k , 



(2vr) 2 - 2e k[ 



a, 



W) 



-ik i b 



1) 



(17) 



The argument of the strong coupling is taken to depend on the transverse momentum k\ in 
order to account for higher-order quantum corrections, originating from momentum scales 
larger than this j2|| |26| . It is obvious that the above integral is not well-defined at very small 
mass scales owing to the singularities of the one (or higher) loop QCD running coupling 
in this region. This makes its evaluation at the edge of phase space sensitive to their 
regularization. The effect of regularizing the k± integration in the infrared region amounts 
to including power corrections to the original perturbative result which are the footprints of 
soft gluon emission at the kinematic boundaries to the non-perturbative QCD regime 

By the same reasoning as applied in the previously considered case, we impose analy- 
tization as a whole, and using Eq. (JSjL we perform the k± integration in (fTTj) to obtain 



[n« (q 2 )1 

with the Borel transform given by 



i 



(18) 



— j sm(7ra/3i) J d£g(£) 



+ 



7T 



b 2 A 2 



where 



9(0 
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d 2 ~ 2e q 1 1 



-2iqb 



(19) 



(20) 



(27r) 2 - 2e q 2 q 2 + £ 

Combining denominators in Eq. fl20l) and carrying out the integrations over we then find 



r n w (q 2 )] = ?l u%\y H doe -^Mvw) 1 

L ;j an n ^ ) J q r(l + rr/3i) 



r (1 - a(3 x - e) 



£ 

n=0 



(n + 1)! V 4 



^2^2 \ n+l-<r/3i-e 



1 

n+1 — a(3\ — e 



:f 



C F Jb 2 k 2 ' 



(21) 



We take the logarithmic derivative because perturbative QCD predicts not the absolute magnitude of 
reaction amplitudes, but only their variation with momentum. The derivative also eliminates the collinear 
divergence related to the integrations over fc ± . 
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where 



'b 2 A 2 



f 



n=0 
1 



n 



si(7r(n + 1)) 



1 ! 



2A2\ n+1 
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6 2 A 2xn+l 



r 



2a2~ 



-n — 1, 



b 2 A 



b 2 A 



n + 1 

oo i 

n=l v 
oo j 



r 



-n -2,1) • 

1 6 2 A 2 
-- m 

2 4 



n 



+ ^(n + l) 



n=l 



1) 



n+1 2si(7rn) A 2 A 



7T 



'.2 A 2 \ " 



1 6 2 A 2 1 1 

--In 1 

2 4 2n+l 



V»(n + 1) 



(22) 



with T(x, y), denoting the incomplete Gamma function [2 



si (x) 



sint 



dt 



2 



x sini 



dt 



(23) 



and ip( x ) — (d/dx) lnT(x). 

At this point some important remarks are in order. The first term in Eq. (|21j) . viz., the 
integral over a, diverges for a Pi = 0, i.e., for small values of a s (k±) (or equivalently for large 
transverse momenta k±). This ultraviolet (UV) divergence is taken care of by the defect of 
the dimension e within the MS renormalization scheme we have adopted. Were it not for 
the terms containing powers of bA, our expression (|21|1 and that found by Korchemsky and 
Sterman in p| (namely, their equation (18)) would be the same. In our case, however, the 
imposition of analytization cures all divergences related to IR renormalons that are generated 
by the T-functions whenever a Pi is an integer different from zero. On the other hand, when 
a Pi is an integer different from zero, say, aPi = m, then the integrand in Eq. (|2ip takes the 
form 



-1) 



m—l 



[ml 



hi 



6 2 A 2 



n=0 



-ir 



in + l)\m\ 



6 2 A 2X«+1 



n 



1 — m 



where 



B = l- lE 



m—l 
+ E k+1 



(24) 



(25) 



k=l 



and the prime on the sum symbol reminds that the term n = m — 1 is excluded. (Obviously, 
for m = 1 the sums above vanish.) 

Let us now turn our attention to the second term in Eq. (|2ip. Evaluating further this 
term and retaining only the leading contributions in 6 2 A 2 , we finally obtain 



h 2 \ 2 h 2 \ 2 h 2 \ 2 
f (b 2 A 2 ) = -a - ai—j- In —j- + + O (6 4 A 4 ) 



(26) 



7 



in which the following abbreviations 



a 







1 si(-7r) 

Z 7T 



3.18 



a 2 = -Ao'Je - ^2 
~ -2.51 



n=0 



+ 



.4 



n+l 



n + 1 (ra + 2)(n + l) 



+ A r(-2,i) + ^(2) + 



2si(V) 



7T 



+ V>(2) 



have been used with 7^ = 0.5772 . . . being Euler's gamma function and 

1 



(n + l)! 



-I)- | 1 | si(7r(n + l)) 



7T 



(27) 



(28) 



Though the integral in Eq. ()21|) (i.e., the first term) cannot be computed in closed form, 
it can be expanded in terms of powers of b 2 A 2 . The only singularity of the integrand is a 
single UV pole at a(3\ = 0, which is regularized dimensionally. For aj3\ an integer, both 
terms inside the bracket have poles, but they mutually cancel so that their sum is singularity 
free and therefore the integral is well defined. Retaining terms of order b 2 A 2 , it is apparent 
from expression (J24j) that the main contribution stems from the leading renormalon a(3\ = 1. 
TIig rGSU.lt is 

[n ,1, (« 2 )] m =[n (I, (0 2 )]p T +[n" , (« 2 )] P o„ (29) 

with the perturbative part being defined by 



[ri (1) (q 2 )] 



C F , In (C/b 2 A 2 ^ 
In 



PT 



TTft In (Q 2 /A 2 



(30) 



where, within the MS scheme, we have set for the renormalization scale (the collinear 
limit) /i 2 = 2Q 2 . Notice that this leading-order result coincides with the one obtained by 
Korchemsky and Sterman (Eq. (20) in 0]). Power corrections in the impact parameter b are 
encoded in the second contribution (6 2 A 2 1): 



[n« (q 2 )] 



pow 



where 



and 



So 



S + b 2 S 2 (b 2 A 2 ) + O (6 4 A 4 ) 



a = 



S 2 {b 2 A 2 ) 



x , , 6 2 A 2 
A ai In a 2 



(31) 
(32) 
(33) 



4tt/V~ 4 

with the constants ao, ai, and a 2 already introduced in Eq. (|2*7j) . 

Hence, the Drell-Yan cross-section V4dy, comprising the leading logarithmic perturba- 
tive contribution (alias the leading Sudakov exponent SW) and including the first power 
correction in 6 2 A 2 reads 



W DY (b, Q) = exp [-S PT (b, Q) - b 2 S 2 (b, Q) + . . .] 



(34) 
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Fig. 1: Exponent of the first non-constant power correction to the Drell-Yan cross-section as a 
function of the impact parameter b. 



S 2 {b, Q) ~ S 2 (b 2 A 2 ) InQ + const 



(35) 



Note that the Q-dependence arises due to collinear interactions, i.e., through the integration 
of Eq. (fTTj) . While the Sudakov factor, representing the perturbative tail of the hadronic 
wave function 0, H3, 2^] , suppresses constituent configurations which involve large impact 
space separations, the exponentiated power corrections in b 2 (leaving aside the constant term 
So), which are of nonperturbative origin, provide enhancement for such configurations, since 
S 2 (5 2 A 2 ) (see Fig. |TJ) is always negative. The consequence is that combining (resummed) 
logarithmic radiative corrections and power-behaved corrections in b, the latter arising from 
soft (nonperturbative) gluon emission and being therefore universal, the net result is less 
suppression of the Drell-Yan cross-section |4| . Here we have an immediate link with the work 
of Akhoury et al. jHf who pointed out within the renormalon approach 0| that the form of the 
exponent in Eq. ()34|) is the same as for the Fourier-transformed pion wave function. In both 
cases the leading power correction in b 2 A 2 has an exponential (Gaussian) form. 2 However, 
and most importantly, with our method we can go beyond their analysis and specify the 
absolute normalization of the power correction that cannot be fixed within perturbation 
theory. In our Drell-Yan calculation, this coefficient, see Eq. (j22I), can be computed explicitly, 
and it turns out to have the opposite sign relative to theirs and depend logarithmically 
on the impact parameter. The upshot is that the inclusion of power corrections leads to 
an enhancement of the pion wave function in b space, counteracting partly this way the 
suppression provided by the familiar Sudakov factor, similar in this respect to the observation 
made in j22J , with the endpoint region bA ~ 1 (where bA is not a small expansion parameter 



2 Let us mention in this context that the Gaussian dependence on the impact parameter b for the Drell-Yan 
process was already noticed by Collins and Soper [30j in their Sudakov analysis. 
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and therefore Eq. ()31|) becomes inaccurate) being less enhanced relative to small b transverse 
distances (cf. Fig. |TJ). 

4. Conclusions 

We have focused on two specific cases - the pion form factor at leading power in Q 2 and 
the Drell-Yan process - which expose all the salient features of the proposed analytization 
methodology. On account of analytization of hadronic functions as a whole, the dispersive 
conjugate of the running coupling is defined unambiguously. Moreover, and even more im- 
portant, one can calculate not only the power of power corrections to hadronic processes, but 
also their concomitant coefficients because this approach does not contain an IR renormalon 
ambiguity from the outset. In this way, we were able to compute explicitly the first power 
correction in to the pion form factor, as well as a Sudakov-type factor to the Drell-Yan 
cross-section which contains the leading power correction in b 2 A 2 . Further applications and 
phenomenological implications of our approach will be pursued in forthcoming publications. 



Acknowledgments 

We wish to thank Alexander Bakulev, Johannes Bliimlein, Christos Ktorides, Sergey 
Mikhailov, Dmitri V. Shirkov, and Wolfram Schroers for stimulating discussions. This 
work was supported in part by travel grants by the COSY Forschungsprojekt Jiilich/Goeke 
(A.I.K.) and Athens University (N.G.S.). 



10 



N. G. Stefanis, Eur. Phys. J.direct C 7 (1999) 1 hep-ph/9911375 . 
D. V. Shirkov, |hep-ph/0012283l 

D. V. Shirkov and I. L. Solovtsov, Theor. Math. Phys. 120 (1999) 1210 [Teor. Mat. Fiz. 120 
(1999) 482] [hep^ph/9909305|. 



G. P. Korchemsky and G. Sterman, Nucl. Phys. B 437 (1995) 415 |hep-ph/94112ll| . 

R. Akhoury, A. Sinkovics, and M. G. Sotiropoulos, Phys. Rev. D 58 (1998) 013011 

|hep-ph/9709497| . 

A. V. Efremov and A. V. Radyushkin, Phys. Lett. 94B (1980) 245. 
G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22 (1980) 2157. 

D. Miiller, Phys. Rev. D 49 (1994) 2525; Phys. Rev. D 51 (1995) 3855 |hep-ph/9411338] . 
S. Agaev, |heiTph/9611215| 



S. Agaev, Mod. Phys. Lett. A 13 (1998) 2637 hep-ph/9805278 ; Eur. Phys. J. C 1 (1998) 321 
|hep-ph/9611283] . 

D. V. Shirkov, Teor. Mat. Fiz. 119 (1999) 55 [Theor. Math. Phys. 119 (1999) 438] 
|hep-th/9810246| ; |hep-ph/0009106 



B. V. Geshkenbein and B. L. Ioffe, JETP Lett. 70 (1999) 161 hep-ph/9906406 . 

A. P. Bakulev, A. V. Radyushkin, and N. G. Stefanis, Phys. Rev. D 62 (2000) 113001 [hep- 

ph/005085]. 

A. V. Radyushkin, JINR preprint E2-82-159, JINR Rapid Communications 4[78]-96, 9 (1982) 
|hep-ph/9907228] . 

N. V. Krasnikov and A. A. Pivovarov, Phys. Lett. 116B (1982) 168. 
A. A. Pivovarov, Nuovo Cim. 105 A (1992) 813. 

D. V. Shirkov and I. L. Solovtsov, JINR Rapid Commun. 2 [76] (1996) 5 hep-ph/9604363 ; 



Phys. Rev. Lett. 79 (1997) 1209 hep-ph9704333 ; Phys. Lett. B 442 (1998) 344 
|hep-ph/971125T] ; D. V. Shirkov, Nucl. Phys. Proc. Suppl. 64 (1998) 106 lhep-ph/9708480] . 
A. I. Karanikas, W. Schroers, and N. G. Stefanis, in preparation. 
J. Qiu and G. Sterman, Nucl. Phys. B 353 (1991) 105; ibid. B 353 (1991) 137. 
H. Contopanagos and G. Sterman, Nucl. Phys. B 400 (1993) 211; ibid. B 419 (1994) 77 
hep-ph/9310313 . 



R. Akhoury and V. I. Zakharov, Phys. Lett. B 357 (1995) 646 hep-ph/9504248 ; Phys. Rev. 
Lett. 76 (1996) 2238 |hep-ph/9512433] . 



M. Beneke and V. M. Braun, Nucl. Phys. B 454 (1995) 253 |hep-ph/ 9506452 . 

G. P. Korchemsky and G. Marchesini, Phys. Lett. B 313 (1993) 433. 

G. Sterman, Nucl. Phys. B 281 (1987) 310. 

G. Curci and M. Greco, Phys. Lett. 92B (1980) 175. 

J. Kodaira and L. Trentadue, Phys. Lett. 112B (1982) 66. 

N. G. Stefanis, W. Schroers, and H.-Ch. Kim, Eur. Phys. J. C 18 (2000) 137 |hep-ph/0005"218] ; 
Phys. Lett. B 449 (1999) 299 |hep-ph/9807298l . 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, 
San Diego, 1980). 

N. G. Stefanis, Mod. Phys. Lett. A 10 (1995) 1419. 

J. C. Collins and D. E. Soper, Nucl. Phys. B 193 (1981) 381; Nucl. Phys. B 197, (1982) 446. 



11 



